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ABSTRACT 

Filamentary structures are ubiquitous from large-scale molecular clouds (few parsecs) to small- 
scale circumstellar envelopes around Class sources (^1000 AU to ~0.1 pc). In particular, recent 
observations with the Herschel Space Observatory emphasize the importance of large-scale filaments 
(few parsecs) and star formation. The small-scale flattened envelopes around Class sources are 
reminiscent of the large-scale filaments. We propose an observationally derived scenario for filamentary 
star formation that describes the evolution of filaments as part of the process for formation of cores 
and circumstellar envelopes. If such a scenario is correct, small-scale filamentary structures (0.1 pc 
in length) with higher densities embedded in starless cores should exist, although to date almost 
all the interferometers have failed to observe such structures. We perform synthetic observations of 
filaments at the prestellar stage by modeling the known Class flattened envelope in L1157 using 
both the Combined Array for Research in Millimeter-wave Astronomy (CARMA) and the Atacama 
Large Millimeter/Submillimeter Array (ALMA). We show that with reasonable estimates for the 
column density through the flattened envelope, the CARMA D-array at 3mm wavelengths is not able 
to detect such filamentary structure, so previous studies would not have detected them. However, 
the substructures may be detected with CARMA D-fE array at 3 mm and CARMA E array at 1 
mm as a result of more appropriate resolution and sensitivity. ALMA is also capable of detecting 
the substructures and showing the structures in detail compared to the CARMA results with its 
unprecedented sensitivity. Such detection will confirm the new proposed paradigm of non-spherical 
star formation. 



1. INTRODUCTION 

It is becoming clear that filamentary structures (few 
parsecs to 10 parsecs in length and typically 0.1 pc in 
width) in molecular clouds are common and need to be 
understood. One clear example is the Integral-Shaped 
Filament region in the north of Orion-A, comprising 
OMC 1-4 , where prestellar cores and p rotostars are form^ 
ing (e.g.,|Chini et al.|1997l |Johnstone fc Bally 19 991 |Aso 
et al. 12000! 'Mutter & Ward- i'iiomps~ ||2007( |Ikeda et al. 

also com ' 
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2007| [Takahashi et al. 2008). Taurus also consists of sev- 
eral large filaments , each with ongoing star-fo rming ac- 
tivity (e.g. Mizuno et al. 1995 Onishi et al. 1998|,Kenyoii| 
et aH2008 1 . I'he large-scale filaments are even more evi- 
dent in recent o bservations from th e Herschel Spa ce Ob 
servaiory (e.g. 



2010| |Arzoumanian et al.||2011[ |Hiir et al.||201f| ), and 



Andre et al. 2010 Men'shchikov et al 



these observations further suggest a tight connection be- 
tween the formation of dense cores and gravitationally 
unstable filaments. While the mechanisms for forming 
these filamentary structures are still under debate (e.g., 

I -., - A ^ ■ .1 . , , ■ I I-. , ' I 



Mac L ow fc Klessen'2004'; 'Heitsc h et al.|2008t |Nakamura| 



& Li|12008; Myers 2009, 2011; Pon et al.||201ip, an ob 



servationally derived processTias been suggested: first 



the filamentary structures at large scales form, possibly 
as a result of magnetic-hydrodynamic (MHD) turbulence 
in the ISM, and secondly, the prestellar cores form from 
the fragments of a subset of filaments through gravita- 
tional instability. 

As the role of large-scale filaments in molecular clouds 
has received significant attention, a number of the lat- 
est observations have unveiled filamentary structure at 
smaller scales (~ few tenths of parsecs in le ngth and 
lOOths of par secs in width). For instance, Hacar & 
Tafalla ( 2011[ ) have observed four subsonic, velocity- 
coherent filanients in L1517 (~0.5 pc in length) that are 
possibly condensed out from the more turbulent natal 
cloud and le ad to the quasi-static fragmentation of cores. 
In addition, Pineda et al. (2011) probed the Barnard- 
5 star-forming core with high angular resolution and 
discovered filamentary structures with ^0.1 parsecs in 
length. The filaments in Barnard-5 are possibly the re- 
sult of fragmentation in a coherent region where subsonic 
motions dominate, and are likely to form stars via future 
gravitational collapse. 

In addition to the filamentary structures in molecular 
clouds at the early stage of star formation, typical length 
of few parsecs and width of 0.1 pc, small-scale filamen- 
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tary structures, typical length of few thousand AU to 0.1 
pc and width of few hundred to few thousand AU, have 
also b een observed in th e envelopes around Class proto- 
stars ( Tobin et al.|2010" ). These filamentary structures in 
the protosteliar envelopes are mostly irregular and non- 
axisymmetric in morphology, suggesting the initial non- 
equilibrium from the prestellar stage. The filamentary 
structure presented near the Class source is reminis- 
cent of the large Herschel observed structures, although 
the size scales of the two are distinct and the properties 
are presumably different. 

The relationship between the large-scale filaments in 
molecular clouds and small filamentary envelopes around 
young protostars still requires further investigation. Sev- 
eral numerical simulations have shown that large-scale 
filaments in molecular clouds a re prone to fragmentation 
leading to p restellar cores (e.g.,|Inutsuka & Miyama 1997 



Hartmann|2002 ) , and filaments a re possibly the most fa- 



vorable mode tor fragmentation (Pon et al. 2011 



2012). 



amen- 



Moreover, studies have also demonstrated that fi 
tary geometries at large scales have a significant impact 
on the geometri es and symmetries of the subsequently 
collapsing cores dSmith et al.]|20TT| ). 

These observations and numerical simulations deliver a 
clear message: filamentary structures from large to small 
scales are clearly playing an important role to the star 
formation process. In this paper, we suggest an obser- 
vational evolution between filaments at the large scale 
and filaments on the small scale. The proposed scenario 
suggests that the small-scale filamentary structure (few 
thousand AU) in protosteliar envelopes originate from 
the filamentary structure (0.1 pc) embedded in the larger 
envelopes of starless cores instead of being produced by 
the protosteliar collapse. We will further show why the 
filamentary structures in starless cores have not been ob- 
served to date. 

2. EVOLUTION OF PROTOSTELLAR STRUCTURE 

It has been well known that dust emission maps 
of Class sources show very spherical emission (e.g., 
Looney et aT][2000l [Shirley et al][2000l [Motte fc Andre" 
20011) . Although molecular surveys of dense cores showed 



non-spherical structures (e.g., Myers et al. 1991), these 
non-symmetric structures were otten considered to be 
material not directly involved in the star formation pro- 
cess, i.e. part of the larger-scale molecular cloud or 
clump, so these components were rarely used in the ob- 
servational modeling of these sources. Instead many au- 
thors assumed that the spherical dust emission indicated 
spherical collapse (e.g., Shu 1977 Terebey et al. 1984) 
and used this symmetry to derive envelope properties 
and place constraints on any e mbedded disk components 



e.g. [K eenc & Masson 1990; Loone y et al. 2003 Har 



vey et al. 2003: J0rgenscn et al. 2009). However, recent 



studies have shown the envelope structures to be more 
complex. 

2.1. Changing the Paradigm for the Inner Envelope of 
Class Protostars 

The ability to use 8 ^m absorption against PAH emis- 
sion background allows the decoupling of the dust den- 
sity a nd temperature for the first time in C lass sources 
(e.g., Looney et al.[|2007 Tobin et al. 




Fig. 1. — Exam ple of the envelope s detected with Spitzer 8 mi- 
cron imaging from |Tobin et alT] | |2010[ | . The images are from Spitzer 
8.0 fim observations. 



of the envelope are complex, filamentary, and often non- 
axisymmetric structures (~ 1000 AU to 0.1 pc). Figure 
[T] illustrates the diversity of structures seen in the Tobin 
sample. IRAS 16253-2429 is what one would expect to 
see in a spherical envelope case. The 8 /im absorption is 
not a good tracer at the central source or in the outfiow 
cavity, since in both cases there is emission in addition 
to the background. In stark contrast, L673 is a clear 



example of the main point of Tobin et al. (2010), which 
is that flattened, filamentary, and non-axisymmetric en- 
velopes are the typical envelope structure. 

How does this result reconcile with interferometric dust 
emission observatio ns which show spher ical emission in 
these sources (e.g., 
tant to remember t 



2000)? It is impor- 



Looney et al. 

lat dust emission depends on both 
dust density and temperature. With flattened or non- 
axisymmetric envelopes and/or outflow cavities in young 
sources, the heating will be inhomogeneous; the lower 
density material near the central source is heated more, 
leading to temperature and density gradients, and the 
dust emission will appear more spherical even if the dust 
distr ibution is not . A go od example, shown in Figure 
from Chiang et al. ( 2010[) , is the source LI 157. Althoug" 
there is a flattened and filamentary envelope detected in 
both N2H+ and the 8 fj,m absorption (also seen in Fig- 
ure the dust emis sion is very spherical and typical of 
a Class protostar. Chiang et al. (2012) constructed a 



measurements, it was realized that 



2010). With these 



he dense portions 



model that has a flattened geometry similar to the N2H 
and 8 /im absorption features and yet still predicts the 
observed spherical dust continuum when non-spherical, 
self-consistent temperature solutions are used. 

However, with enough sensitivity the filamentary 
structures can still be seen in dust emission. Figure |3] 
is the dust emission toward L1157 with the Submillime- 
ter Array (SMA) at A = 1.3 mm (Tobin et al. 2012, in 
prep). In this case, they detected the extension along 
the flattened envelope and even an extension along the 
outflow (also see Stephens et al. 2012, in prep). The 
extension along the outflow illustrates how the heating 
is facilitated by lower density material (in this case in 
the outflow cavity). In other words, the heating in these 
sources are not uniform, which can lead to a distortion 
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Fig . 2. — Example of dust and N2H+ differences from [Chiang| 
|et al.| ( ["2010[ l. However, both emission was fit in a simple density 
model. 
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Fig. 3. — With higher sensitivity, SMA observations detect the 
extension of the flattened envelope and the increased heating along 
the outflow cavity. The blue and red arrows indicate the outflow 
directions of the blue-shifted and red-shifted component, respec- 
tively. The beam size is in the bottom-right corner. The contours 
are ±2, ±3, 6, 9, 12, 22, 30, SOct, where a = 2.16 mjy beam-i. 



in the structure suggested by only the dust continuum. 

Indeed, when comparing the observations of L673 and 
L1157 with the traditional view of spherical star forma- 
tion, we need to change the cartoons of star formation. 
Figure |4] demonstrates our suggestion of moving from 
spherical star formation structures to filamentary star 
formation structures in Class protostars to be more 
consistent with observations. The left panel presents the 
traditional model assuming sphericity that has impacted 
our theoretical understanding for decades. In this model, 
protostellar collapse is axisymmet ric and sph erical based 
on a singular isother mal sphere (|Shu[ 1977), With the 
inclusion of rotation ( |Terebey et aT 19M), the density 



structure is slightly flattened and mostly remains spher- 
ical beyond the centrifugal radius. On the other hand, 
the right panel in Figure [4] shows the axisymmetric and 
filamentary envelo pes that are often seen in our Class 
observations (e.g., Tobin et al. 2010). The filamentary 
envelopes with higher density are forming inside the am- 
bient cloud at lower density. 

The change from spherical view of star formation to 
filamentary view certainly has important consequences, 
as several analysis techniques are based on the assump- 
tion of sphericity. For example, for a single beam mea- 
surement of the region the spherical assumption will sig- 
nificantly underestimate the mean density of the dense 
material. Alternatively, low sensitivity interferometric 
maps will concentrate on the peak and likely miss the 
large structure and thus the shape. In addition, blue- 
skewed spectra have been extensively observed with op- 
tically thick molecular lines in starless cores. The in- 
terpretation of the blue asymmetry, together with opti- 
cally thin lines peaking in the absorption dip, has been 
spherical collapse. Moreover, the modeling of spectral 



energy dis tributions ( SEDs) (e.g., [Whitney et al 2003 



[Robitaille et al. 2006) extensively used m Class 6 anc 



Class 1 sources is based on spherical/ axisymmetric mod- 
els. As non-spherical envelopes are more common, spher- 
ical models may not provide accurate descriptions of pro- 
tostellar properties, so extra caution needs to be applied. 

2.2. Observationally- Driven Scenario for Filamentary 

Collapse 

Filamentary structures appear to be ubiquitous from 
large molecular clouds to small scale circumstellar en- 
velopes. These filamentary structures are also observed 
to be tightly connected to the star formation process, as 
prestellar cores and young protostars are located within 
these filaments. From these observations, we propose an 
observationally derived scenario of filamentary collapse 
in star formation that is summarized in Figure [5] As 
shown in the cartoon, there are approximately five steps 
in our observational-based picture of the filamentary col- 
lapse process. Among the five steps. Step I, II, IV and 
V are from observations, and Step III is a prediction of 
high-density filamentary structures in starless cores, to 
connect Step II and IV. 

In Step I, molecular clouds are formed as filaments 
with a few parsecs to 10 parsccs in leng th and a char 
acter istic width of 0.1 pc (S chneider & Elmegreen||1979 
Ball y et al.]|T987| [Johnstone fc Baliy 1999; Andre et a 
|2010| " These large-scale filaments are probably turbu 
lent and prone to fragmentation, leading to subsequent 
velocity-coherent, higher density filaments (few tenths of 
parsecs in length) that are considered as the birthplaces 
of pre stellar cores ( Inutsuka fc Miyama[[1997| [Hartmaim 



ir- 

3] 



2002 Pon et al.|20TT 



T] ). In Step 11, the fragmented fila- 
ments collapse along tne long and short axis while feedin g 



material along the filament (e.g., Hacar fc Tafalla|2011 ) 
enhancing the mass in a location and forming a higher 
density oblate (or prolate) starless core as observed with 
single-dish o bservations (,Curry^2002j_ Jones fc Basu^2002[ 



Tassis 2007 Tassis et al.|2009). I'he core tormation. may 



be related to the hows from large-scale motions along the 
larger filaments and is kinematically coupled with the 
parental cloud ( Hacar fc Tafalla[[2011 ). This would form 
a higher density filamentary structure embedded inside of 
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Fig. 4. — Schematic view showing the changing view of star formation from spherical collapse (left panel) to filamentary collapse (right 
panel). Note that the size scales are exaggerated to better illustrate the structures. 



the starless core as seen in Step III. This substructure is 
the kinematic descendant of the flow along the larger fil- 
ament and the origin of the filamentary envelopes seen in 
the Class objects. As the collapse continue s in Step IV, 



mate rial infalls along the smaller filament (Tobin et al 
2012 1 and the oblate (or prolate) starless core contmues 



to collapse into a centrally condensed envelope of a Class 
protostar. A Class source is created (~ 5000 AU in 
size), while the large scale filamentary structure (~ 1000 
AU to 0.1 pc) remains behind containing an appreciable 
fraction of the total mass of the envelope plus source. In 
Step V, the Class source evolves to a Class I source with 
a protostcllar disk and the larger structure dissipates. 

Although this scenario fits together, there is one se- 
rious problem with our proposed evolution of filaments 
in star formation: no one has detected the substructure 
(e.g. filamentary structure) predicted in Step III in star- 
less cores to date. Unfortunately, there is some difficulty 
in detecting these structures. One could use molecular 
line tracers such as N2H"'" or NII3, which often corre- 



spond to the 8 /i m absorption (e.g., Chiang et al.||20l0 
Tobin et al. 201 1 1 . However in starless cores, 1N2H "'" couk _ 



still ha ve chemical effects such as depletion (e.g., iBergin 
et al.|2002), although several studies showed I'l'afalla et al 



(200ii) less depletion for N2II+ than other molecules 
Since the depletion usually ocurrs at the center of the 
core, the filament could appear fragmented in the map 
depending on the size of the central depletion. In ad- 
dition, the molecular distribution could originate from 
chemistry and not well-trace the dense material. Thus, 
to confirm detection of substructure, we must rely on 
dust continuum emission. Dust emission at millimeter 
wavelengths presumably is more appropriate than 8 /im 
extinction because of the low optical depth. 8 /im ex- 
tinction shows detection in the outer regions only if the 
background signal to noise is high enough. 

In order to resolve the structures, we must have reso- 
lution of ~ 5 arcsec, which im plies interferometers. For 



observations at 3 mm toward 11 starless cores in Perseus 
with CARMA. Although two sources were d etected, they 
were later reclassified as protostcllar objects (Enoch et al. 
2010 Schnee et al. 2012 ), implying only non-detections ot 
sub-structure of starless cores, contrary to our suggested 
evolutionary sequence. Our explanation is that the sub- 
structure was not detected due to a lack of sensitivity. In 
the following section, we investigate that possibility and 
place con straints on the unde rlying filamentary structure 
based on Schnee et al. (2010) results. 



example, Schnee et al. (2010) performed dust continuum 



3. SYNTHETIC OBSERVATIONS 

To examine the likelihood of our proposed struc- 
tures in starless cores, we make synthetic observations 
with CARMA, dire ctly comparing to the observations of 
Schnee et al. (2010), and ALMA, using the fiattened en- 
velope around L1157 as a model. The result will show 
that the expected structures are below CARMA D ar- 
ray's detection threshold at 3mm, but they should be 
detectable with CARMA D and E array observations at 
3mm, CARMA E array at 1mm, and ALMA 1mm obser- 
vations. This implies that there is not yet a disagreement 
between our observational-based proposed evolutionary 
scheme in Step III for low-mass star formation and cur- 
rent observations, and an exciting observational future is 
suggested. 

3.1. CARMA observations 

We simulate CARMA imaging with parameters used 
by Schnee et al. (2010): heterogeneous array (six IO- 
meter antennas and nine 6-meter antennas) imaging with 
the CARMA-D array configuration at 3 mm continuum. 
We use the Miriad tasks uvgen, demos and uvmodel, 
based on |Wright (2010) without the 3.5-meter telescopes. 
To find the detection limit of such a structure, we also 
simulate CARMA D-|-E array observations at 3 mm and 
CARMA E array at 1 mm. Baselines range from 3fcA to 
38fcA for the D array at 3 mm, 2fcA to 19A:A for the E 
array at 3 mm, and 5fcA to 47fcA for the E array at 1 
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Fig. 5. — Illustration of filamentary collapse in five Steps. In 
Step I, molecular clouds are formed with filamentary shapes (few 
parsecs in length and 0.1 pc in width) that are prone to fragmenta- 
tion. In Step II, the subsequently-fragmented filaments (few tenth 
of parsecs in length) collapse and form prolate or oblate starless 
cores (0.1 pc in size). In Step III, embedded in the starless core is 
a filamentary structure of higher density that arises from the fiow 
along the the large-scale filament axis. In Step IV, the starless core 
continues to infall into a centrally condensed envelope of a Class 
protostar (~ 5000 AU). In Step V, the Class source evolves to 
a Class I source with a protostellar disk (few hundred AUs) and 
outflow (the arrows). The orientation of the protostellar disk de- 
pends on the detailed kinematics of collapse and is not necessarily 
along the fllament as shown in this cartoon. Note in each step the 
structures with lower density are indicated with dashed lines. 



mm. The observing rest frequency is centered at 90 GHz 
for 3 mm observations and 230 GHz for 1 mm observa- 
tions with a total bandwidth of 4 GHz for continuum 
observations. The total observing time on the target is 
6 hours for each synthetic observation (for the one with 
CARMA D+E array, the observing time is three hours 
for the D arr ay an d three hours for the E array). In the 
analysis (Sect 3.31, we perform a small mosaic (standard 



seven-pointings) around the source to capture all the ex- 
tended structure. Table [l] summarizes the synthesized 
beam sizes and noise levels for the simulated CARMA 
observations. 



3.2. ALMA observations 



We used the task sim_observe and sim_analyze in the 
package casapy to perform the simulated observations 
with ALMA. The angular resolution is requested to be 
1.2" in the simulation to observe detailed structures, and 
a small mosaic is applied to capture all possible structure. 
The observing time for each mosaic pointing is 100 s, 
and the total observing time is 2 hours. The observing 
frequency is centered at 90 GHz with the bandwidth of 8 
GHz for continuum observations. Thermal noise is added 
with a typical precipitable water vapor of 2.8 mm. The 
clean threshold is set to 1.5 times the noise rms, and the 
pixel size is set to 0.12 arcsecs. 

3.3. Modeling and Results 

As decribed in Step HI and IV in Figure[5] the flattened 
protostellar envelopes around Class sources are specu- 
lated to be highly connected to the filamentary structure 
at the previous stage, the prestellar phase. Therefore, 
the envelopes around Class sources best describe the 
morphology of the high-density filamentary structures in 
Step HI. To simulate the structure in Step HI, we mod- 
ify the flattened envelope in the Class source LI 157 by 
the physical conditions expected at the prestellar stage, 
in order to examine if CARMA and ALMA are able to de- 
tect filamentary structure at the prestellar stage. Figure 
pi a) shows the ext inction ma p of LI 157 from the Spitzer 
8/im observation (Looney et al. 2007). We chose L1157 
to model as it has an obvious tilanicntary envelope struc- 
ture seen in the 8/Ltm absorption against the background 
emission, and the symmetric structure c an be approx- 



imate d with a radial density power-law (Looney et al. 
20071. The spatial scale of the filamentary envelope is 



U.l p c, too large for a c ircumstellar disk or a pseudo disk 
(e.g., |Gani fc Shu|1993[ ). The L1157 dark cloud is located 
^ 250 parsec^away with an edge-on view concealing the 
Class source embedded in the flattened envelope nearly 
perpendicular to a large powerful outflow from the north 
to the south. The distance to L1157 is approximately the 
same as Perseus (also at ~ 250 pc) and thus rem ains an 
excellent pro xy for Perseus when comparing with |Schnee| 
et al. (2010). If the proposed scenario of filamentary 
collapse IS correct, the flattened envelope is expected to 
be related to the filamentary structures on larger scales, 
and thus this source is suitable for modeling the transient 
phase in the prestellar stages with appropriate physical 
conditions. 

To better concentrate on the filamentary envelope, we 
removed the emission from the outflow and the scattered 
light from the central object. We then filled the inner 
regions with the averaged value from the envelopes on 
the two sides, as shown in Figure lejb). The total mass 
calculated from the extinction increased by about 10% by 
filling the inner region with this method. The extinction 
map was compared toward background stars measured 
in the near-IR to an optical-depth i mage generated w ith 
or without the zodiacal correction ( Tobin et al.||2010 |. 

We next generated the brightness map at millimeter 
wavelengths (at 3 mm in our model) by assuming that 

^ The LI 157 cloud is estimated to have a similar galactic lat itude 
as th e absorbing clouds with 200 pc and 300 pc in Cepheus ( |Kun| 
|1998| and ther efore we adopted a distance of 250 pc in this paper, 
in comparison, |Kirk et al.| j2009| adopted a distance of 325 pc for 
the region around L,ii5Y. 
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TABLE 1 

Synthesized beam sizes and noise levels for. the synthetic observations 



CARMA D array 


CARMA D+E array 


CARMA E array 


3 mm 


3 mm 


1 mm 


L1157 modeling (Scct[33| 5.14" X 4.76" 


7.92" X 7.23" 


4.21" X 3.65" 


0.3 mjy beam^^ 


0.15 mJy beam~^ 


0.35 mJy beam~^ 




1 1 TW^SWiTTiITi ' 




20^39^15' 12' 



9= 6' 3=^ 
a (J2000) 



20"39"'24" 18' 12" 6' 1,9'"0' 5+' +8" 
a (J2000) 

Fig. 6. — (a): the left panel. It is the extinction map of L1157 from the Spitzcr 8 fim observation. A flattened envelope is seen in 
absorption against the background emission around the central Class object with the outflow nearly perpendicular to it. (b): the right 
panel. The emission from the outflow and the scattered light from the central object are removed. The central part is filled with the 
average of the envelopes on the two sides. The color scale shows optical depth. 



LI 157 is optically thin at 3 mm. We first calculated the 
mass contained in each pixel from the extinction map, 



M = dVtx 



1.496 X 10^^ X X 
AU/ 



where d^l is the pixel solid angle (1.2")^, D is the distance 
in parsecs (250 pc for LI 157) and Kg^m the dust plus 
gas opacity at 8^m. The 3 mm flux in each pixel is then 
determined from the mass, opacity, and temperature. 



F = 



M X B^{T) X K3r, 
D2 



where B^, (T) is the Planck function and Ksmm is the dust 
plus gas opacity at 3mm. 

We used 0.00169 cm^ f or K^mm by assumin g 100 
for the gas to dust ratio (e.g., Schnee et al. 2010). The 
temperature was a ssumed to bea con stant at 10 K for 
starless cores (e.g., Schnee et al. 2009[ ). After obtaining 
the brightness map at 3 mm, we simulated the CARMA 
D array observations with our model. Since the opac- 
ity in the infrared is poorly constrained, we generated 
models with varying values for Kspm, which also modi- 
fies the derived mass in the flattened envelope structure. 
The equation above indicates that the observed millime- 
ter brightness decreases with increasing Ksfj.m, since less 
mass is required to produce the IR extinction. In the left 
panel of Figure [7| we show that with an expected value 



of ^ 10.96 c m^ g~i for Kg^^ ( |Tobin et al^|2010[ [Butler 
fc Tan||2009[ ) toward L1157, the CARMA-D array is not 



able to detect the filamentary structures at the prestellar 
stage. For the structures to be clearly deleted (the right 
panel of Figure [t]) , the dust opacity at 8 fim would have 
to be an unphysically small value. 

To fully explore CARMA's capability, we performed 
the synthetic observation with CARMA D-l-E array at 
3 mm since the E array is more compact and sensitive 
to emission at larger scales than the D array. The total 
observing time is six hours (three hours with the D array 
and three hours with the E array). The value for Kg^m 
used is 10.96 cm^ g~^ to compare with the result from 
the D array only, since it produces the weakest emission 
(contains the least mass). As shown in Figure [s] the fil- 
amentary structure is detected with a similar noise level 
(0.28 mJy/beam) as the D array (0.3 mJy/beam), al- 
though only the structures with stronger emissions close 
to the center could be seen and the structures are not 
in detail. The detection suggests that the non-detection 
with the D array is due to a combination of spatial res- 
olution and sensitivity. Furthermore, we shift the ob- 
servation from 3 mm to 1 mm a ssuming that the dust 
opaci ty is 0.9 cm^ g~^ at 1 mm (Ossenkopf & Henning 
1994[ ) and KgMm is still 10.96 cm^ g-\ Again, CARMA E 
array is able to show detection on the structure as shown 
in Figure |9j as the brightness increases toward the short 
wavelengths. 
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Fig. 7. — Simulated observations with CARMA D array at 3 mm for six hours with different values of ksm'"- The center is positioned 
at 20''39'"05.2* (RA) and 68°02'15.3" (Dec) (J2000). The synthesized beam size is 5.14"by 4.76 "shown on the bottom-right corner. The 
noise level cr is 0.3 mjy/beam, and the contours are ±3,±4.2,±6, ±8.5, ±12, ±17, ±24, ±34 x a (in step of y^a). The color scale shows 
flux in Jy/beam. With the reasonable value for Kg^m (10.96 cm-^ g""*^) in the left panel, no structures are detected. For the structures to 
be clearly detected (the right panel), ftg^m needs to be an almost impossibly small value (2.0 cm^ g~^)- 
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Fig. 8. — Simulated observations with CARMA D+E array at 3 
mm for six hours in total (three hours for the D array and three 
hours for the E array). The value for fcgfim is 10.96 cm^ g~^. The 
center is positioned at 20''39™05.2'' (RA) and 68°02'15.3" (Dec) 
(J2000). The synthesized beam size is 7.92"by 7.23"shown on the 
bottom-right corner. The noise level a is 0.28 mJy/beam, and the 
contours are ±3, ±4, ±5, ±6, ±7, ±8, ±9, ±10 X cr. The color scale 
shows flux in Jy/beam. 



We further ran the synthetic observations with ALMA 
also for the case of ks^™ = 10.96 cm^ g~^. Figure 
10 shows the result from the synthetic observation (the 
contours are plotted by percentages of the peak flux in- 
stead of the noise levels due to the artificial effects from 
resolving-out large structures). As can be clearly seen, 
ALMA is able to detect most of the structure in the flat- 
tened envelope. Comparing the results from CARMA 
and ALMA for Kg^™ = 10.96 cm^ g"\ ALMA's un- 
precedented sensitivity greatly improves the appearance 
of filamentary structures and should provide a powerful 
tool for uncovering any hidden filamentary profiles at the 
starless/prestellar stage. Located in the southern hemi- 
sphere, ALMA is not able to look at L1157; these results 




20 
offset (") 



Fig. 9. — Simulated observations with CARMA E array at 1 mm 
for six hours in total. The value for Kg^m is 10.96 cm-^ g""*^- The 
center is positioned at 20'"39™05.2'' (RA) and 68°02'15.3" (Dec) 
(J2000). The synthesized beam size is 4.21"by 3.65"shown on the 
bottom-right corner. The noise level a is 0.49 mJy/beam, and the 
contours are ±3, ±4, ±5, ±6, ±7, ±8, ±9, ±10 X a. The color scale 
shows flux in Jy/beam. 



are, however, indicative of the ALMA observations with 
other starless cores. 

4. CONCLUSION 

In this paper, we posit an observationally derived sce- 
nario for filament-driven star formation that incorporates 
the evolution of star-forming cores with filaments into fil- 
amentary envelopes from large to small scales. Molecular 
clouds are formed as filaments (few parsecs to 10 parsecs) 
and then fragment to smaller filaments (few tenths of par- 
secs), which eventually collapse to form triaxial starless 
cores. As collapse continues the material infalls along the 
filament into a centrally condensed filamentary envelope 
of a spherical Class source, which keeps evolving to a 
Class I source with a protoplanetary disk. 
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cores) and thus we will be able to conduct quick and ef- 
ficient surveys of the geometry of the envelopes around 
starless cores. The scenario proposed scenario can be 
immediately tested by observations with the current in- 
struments. 



-40 -20 20 40 

X offset (") 

Fig. 10. — Simulated observations with ALMA for KS/jm = 10.96 
cm^ g~^. The noise level a is 0.15 mjy beam~^. The contours 
indicate ±10, ±15, ±20, ±25, ±30, ±40, ±50, ±60 X cr . The color 
scale shows flux in Jy/beam. The synthesized beam size is 1.2" 
shown on the bottom-right corner. 

If such a scenario is correct, the filamentary structures 
at the prestellar stage should exist. The only reason 
that they have not yet been detected is sensitivity to 
large-scale emission in the surveys. They are possible to 
detect with CARMA D+E array at 3 mm due to more 
appropriate resolution and CARMA E array at 1 mm as 
a result from the higher brightness at 1 mm; however, 
ALMA is even more capable of clearly detecting detailed 
structures of the filamentary envelopes. In fact, the very 
high sensitivity of ALMA will allow for much shorter in- 
tegrations (less than two hours for L1157-like prestellar 
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